SPECIFIC HEAT OF Y^Lui-^BsC IN THE MIXED STATE 
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The temperature and magnetic field dependences of the electronic specific heat c v in the supercon- 
ducting mixed state as well as the upper critical field H C 2(T) have been measured on a series of 
polycrystalline YzLui-^IS^EhC samples with the aim to study the influence of disorder. The elec- 
tronic specific heat contribution c esm = ^y(H) ■ T in the mixed (vortex) state exhibits significant 
deviations from the usual j(H) oc i/-law for conventional superconductors for all compositions 
resulting in a disorder dependent negative curvature of 7 oc H 1-13 ^. For the samples with sub- 
stitutional disorder at the rare earth site f3 and the positive curvature exponent a of H C 2(T) are 
reduced but a, /3 > holds for all x, as distinct from the case of much stronger disorder due to iso- 
electronic partial Pt substitutions at the Ni site where a and f3 rapidly vanish (Nohara et al. (1999) 
J. Phys. Soc. Jpn. 68 1078). a exhibits a clear minimum near x — 0.5. The magnitudes of H^{T) 
and T c , the curvatures of H C 2(T) and 7(-ff), as well as the Sommerfeld constant 7at in the normal 
state are reduced in a similar fashion for intermediate composition. 



1. Introduction 

The recent discovery of the rare earth transition metal borocarbide (nitride) family 
(RC)„T 2 B 2 C(N), R=Y, Lu, Sc, Th, La; T=Ni, Pd, Pt; n=l to 3 [1, 2] which con- 
tains superconductors with relatively high transition temperatures T c up to 23 K 
has stimulated numerous studies of their thermodynamic and transport properties 
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in the superconducting as well as in the normal state. At first glance most of those 
results support a classification of these materials as typical intermetallic phonon 
mediated superconductors with a moderately strong coupling strength. However, 
clean RM2B2C samples exhibit also some features unexpected for ordinary super- 
conductors both in the superconducting and in the normal state. Among several 
peculiarities [3] we emphasize the unusual shape and the strong disorder depen- 
dence of the upper critical magnetic field H C 2 (T) and the T 3 dependence of the 
electronic specific heat c es (T) in the superconducting state compared with ex- 
ponential behaviour for ordinary s-wave superconductors. According to Nohara 
et al [4] isoelectronic transition metal substitution affects strongly the field de- 
pendence of the electronic specific heat contribution c esm = 7(H) ■ T in the 
mixed state. While for an Y(Nio.sPto.2)2B2C single crystal j(H) oc H, a marked 
sublinear law was observed for pure YM2B2C single crystal and polycrystalline 
LuNi 2 B 2 C[5] 

j(H)/j n oc 

where 7^ denotes the Sommerfeld constant in the normal state. Thus already 
small substitutions in the TB network (which is commonly believed to be the 
crucial locus for the superconductivity) affect significantly the superconducting 
properties, pressumably due to the predominant role of the T-derived d states in the 
mechanism of superconductivity. At the same time by isoelectronic substitutions 
in the RC charge reservoir much weaker disorder could be studied (see Fuchs et 
al, Drechsler et al, these Proa). The remarkable j(H) oc v^-law for YM2B2C 
and LUM2B2C was regarded at first as evidence for (i-wave pairing [5]-[7]. In this 
context it is interesting to study also the closely related Y x Lui_ x Ni2B2C system 
as an appropriate weak disorder case. Thus changing the composition x, more 
detailed information should be gained about the influence of weak disorder upon 
the field dependence of the specific heat c p (T, H), the shape and the magnitude 
of H C 2(T), as well as about the nature of the pairing state. 

2. Experimental Details 

Polycrystalline Y. E Lui_ ;E Ni2B2C samples were prepared by a standard arc melting 
technique. Powders of the elements were weighted in the stoichiometric compo- 
sitions with a surplus of 10 wt.% boron to compensate the high losses of boron 
caused by the arc melting. The powder was pressed to pellets which were melted 
under argon gas on a water-cooled copper plate in an arc furnace. To get homo- 
geneous samples, they were turned over and melted again four times. After the 
melting procedure the solidified samples were homogenised at 1100°C for ten 
days. The specific heat was measured in the range 4.2K < T < 20K and for 
magnetic fields H < 8T using a quasi adiabatic step heating technique. The heat- 
ing pulses were generated by a strain gauge heater and T was measured with an 
Au-Ge thin film resistor. 
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Figure 1. Zero magnetic field specific heat c p (T) /T vs. T 2 of a Y^Lui-^r^EhC series. 
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Figure 2. Specific heat c p (T, H) /T vs. T 2 of YNi2B2C for various magnetic fields. 



3. Results and Discussion 

The specific heat c p /T vs. T 2 at H = of a Y^Lui-ziN^E^C series is shown in 
Fig. 1 . The corresponding curves for H < 8T of the pure Y sample (x = 1) are 
depicted in Fig. 2. Applying H = 8T resulted in a shift of T c below the lowest 
measured temperature of 4.2K. Measurements in this field were used to analyze 
the normal state specific heat c p = jnT + PdT 3 , where 7jv is the Sommerfeld 
constant and /3dT 3 is the Debye contribution. The 7^ was determined by extrapo- 
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Figure 3. Magnetic field dependence of the specific heat contribution j(H) of the vortex core elec- 
trons in the mixed state (H < H C 2) in units of the effective upper critical field H C 2,eff(0) shown 
in Fig. 4. The straight reference line corresponds to the usual s-wave dirty limit behaviour. The 
inset shows the curvature parameter (3 of -y(H) according to 7(iJ)/7jv oc (H/ H C 2, e f /(0)) 1_/3 . 
Noteworthy, the case (3 = 0.31 to 0.35 and H/H c2 , e ff(0) < 0.75 is not very distinct from the dirty 
(unitary) d-wave limit where f(H) oc H In H according to Ref. [17] (not shown here). 



lating the c p /T vs T 2 curves to T = in the normal state. This way we ob- 
tained 7A r= 20.4 (x = 0), 19 (x = 0.25), 18.3 (x = 0.5), 18 (x = 0.75) and 
20.2mJ/molK 2 (x = 1) for our Y-rLui-^rN^E^C series in good agreement with 
the data reported by several groups [9]-[13]. To determine j(H), the c p /T vs. T 2 
curves have been extrapolated linearly to T = from the data in the range 4.2K 
< T < 7K or up to the onset of the transition to the normal state. The results are 
shown in Fig. 3. Apparently, for all Y^Lui-zI^E^C samples j(H) is a sublinear 
function on H. Generalizing Eq. (1), the data were analyzed by the expression 

7 (H)/ 7N = (H/H^eff^x)) 1 -^, (2) 

where the fitting parameter f3 measures the negative curvature of ^y(H) while 
#c2, e //(0) « (0.8 to 0.9)F c2 (0) yields a lower bound for the true H c2 {0) (see 
Fig.4). We obtained p = 0.66, 0.42, 0.42, 0.41 and 0.56 from x = to x = 1, 
respectively. The behaviour of (5{x) is shown in the inset of Fig. 3. The curva- 
ture of 7(.ff) becomes maximal for the bordering cases x = 0; 1 and markedly 
smaller in between. Noteworthy, our curvatures for LUM2B2C and YM2B2C do 
slightly exceed those derived from Eq. (1) as reported in Refs. [4, 5]. Anyhow, 
j(H) shows for all Y. E Lui_. E Ni2B2C samples a clear H sublinear dependence. 
A j(H) w \/^-law as well as the c es (T) oc T 3 dependence raise the question 
whether an unconventional pairing mechanism is responsible for this peculiarity. 
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According to Ref. [6] j(H) oc \fH is a signature for a nodal order parameter 
with (i-wave pairing (for which (3 = 0.57 has been predicted in Ref. [14]) while 
j(H) oc H points to nodeless s-wave pairing in the dirty limit since there are 
conventional superconductors which exhibit also deviations from the j(H) oc H- 
law in the clean limit, e.g. V 3 Si [15], NbSe 2 [4] {(3 = 0.33), and CeRu 2 [14]. 
Hence, the sublinear j(H) dependences observed here are not inevitably related 
to d-wave pairing [16]. 

Since the lattice constants vary linearly between x = and x = 1, one 
might expect a maximum of T c around x = 0.6 according to an "universal" curve 
T c = T c (d,Ni-Ni) proposed by Lai et al. [18], where d^i-Ni is the Ni-Ni distance. 
However, in accord with Refs. [8,9] a dip near x = 0.7 is observed (T c rs 14. 6K 
at x = 0.75). A similar behaviour is found for H C 2(0) and also for the curvature 
exponent a (see Fig. 4): 

H c2 (T) = H* c2 (l-T/T c ) 1+a , valid for 0.3 < T/T c , (3) 

(distinct from a statement of a constant curvature [9]). H* 2 ~ 1.1 -r l.2H C 2(0) 
gives an upper bound for the true value H C 2(0) at T = 0. Ongoing from pure 
samples in the clean limit to the mixed specimens which remain still in the quasi- 
clean limit, all these quantities are somewhat reduced due to weak disorder at 




0.1 b — ■ — 1 — ■ — 1 — ■ — 1 — ■ — 1 — ■ — u 

0.0 0.2 0.4 0.6 0.8 1.0 
Y-concentration x 



Figure 4. Composition dependence of the transition temperature T c derived from the onset of the 
jump of cp(T) (upper panel), the Sommerfeld constant 7jv (second panel), the upper (lower) bound 
for the upper critical field (third panel, see Eqs. (2,3) and text), the specific heat curvature exponent 
P of -y(H) (fourth panel), and the curvature exponent a of the upper critical field H C 2 (lower panel). 
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the R-sites as well as due to a slight reduction of the electron-phonon coupling 
constant A at intermediate x deduced from the similar behavior of 7^ [19, 20] 
(see Fuchs et al, these Proa). 

To summarize, we have shown that the j(H) curves of the pure specimens 
(x = 0; 1) of Y-eLui-zI^E^C exhibit the highest negative curvatures. Weak 
disorder effects, caused by isoelectronic substitutions of Y by Lu, yield a decrease 
of the curvature without reaching the linear limit, which would be the hallmark 
for the dirty limit. Similar moderate suppressions which are typical for the quasi- 
clean limit have been found for the upper critical field, its curvature, and T c . 
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The temperature and magnetic field dependences of the electronic specific heat in the superconduct- 
ing mixed state as well as the upper critical field have been measured on a (Y,Lu)Ni2B2C series to 
study the influence of disorder. The electronic specific heat in the vortex state shows deviations 
from the usual gamma (H) proportional to H-law for conventional superconductors for all samples 
resulting in a disorder dependent negative curvature of gamma proportional to H sup(l-beta). For 
samples with disorder at the R site beta and the positive curvature exponent alpha of H sub(c2) 
are reduced but alpha, beta remain always positive, as distinct from much stronger disorder due to 
Pt substitutions at the Ni site where alpha, beta rapidly tend to 0. Alpha shows a minimum near 
x= 0.5. H sub(c2)(0), Tc, the curvatures of H sub(c2) and gamma(H), as well as the normal state 
Sommerfeld constant gamma sub(N) are reduced similarly. 

Keywords: Specific heat, upper critical field, dirty limit, clean limit, Sommerfeld constant 



